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A solution to the subdiffusion-efficiency paradox: Inactive states 
enhance reaction efficiency at subdiffusion conditions in living cells 
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Abstract -Macromolecular crowding in living biological cells effects subdiffusion of larger bio- 
molecules such as proteins and enzymes. Mimicking this subdiffusion in terms of random walks on 
a critical percolation cluster, we here present a case study of EcoRV restriction enzymes involved 
in vital cellular defence. We show that due to its so far elusive propensity to an inactive state 
the enzyme avoids non-specific binding and remains well-distributed in the bulk cytoplasm of the 
cell. Despite the reduced volume exploration capability of subdiffusion processes, this mechanism 
guarantees a high efficiency of the enzyme. By variation of the non-specific binding constant 
and the bond occupation probability on the percolation network, we demonstrate that reduced 
non-specific binding are beneficial for efficient subdiffusive enzyme activity even in relatively small 
bacteria cells. Our results corroborate a more local picture of cellular regulation. 



Introduction. — Diffusion-limited biochemical cel- 
lular reactions underlying signalling and regulation pro- 
cesses have traditionally been investigated at dilute 
solvent conditions PQ. The relevance of this picture for 
diffusion control in living cells has been challenged in view 
of macromolecular crowding, the occupation of a consider- 
able volume fraction / of the cellular cytoplasm by larger 
biopolymers [2]. Estimates for / typically range from 35% 
to 40%. Bearing in mind that on a cubic lattice the site 
percolation threshold is / ps 31% and that of bond percol- 
ation / w 25% (4j[5j, molecular crowding may indeed ap- 
pear severe. Crowding effects changes in enzyme function 
and turnover, as well as protein folding and aggregation 

mm- 

Larger biopolymers and tracers in living biological cells 
and artificially crowded control environments perform sub- 
diffusion of the form I7hT51 



(r 2 (f)} ~ t a with < a < 1, (1) 
as observed experimentally for particles as small as 10 kD, 
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with a in the range of 0.40 to 0.90 0I3QI1QMI] , in ac ~ 
cord with recent high-detail simulations [5] . The observed 
subdiffusion has been measured to persist over tens to hun- 
dreds of seconds [9l[TT| and thus appears relevant to cellu- 
lar processes such as gene regulation or molecular defence 
mechanisms. Subdiffusion leads to reduced global volume 
exploration, to dynamic localisation at reactive interfaces 
[18j . and may prevent chromosomal mixing in cukaryotic 
nuclei [10] . It has been argued that subdiffusion may in 
fact be beneficial for cellular processes, by increasing the 
probability to find, and react with a nearby target [9"lll6j. 
In that sense subdiffusion would give rise to a more local 
picture of diffusion-limited biochemical reactions in biolo- 
gical cells, see below. 

Here we study by simulations the reaction dynamics 
of EcoRV restriction enzymes in Escherichia coli (E.coli) 
bacteria. EcoRV in solution forms homodimcrs of molecu- 
lar weight 58 kD [53] , and thus belongs to the range of sizes 
for which subdiffusion ((TJ) under crowding was reported 
6,8,131 Our results show that despite the subdiffusion 
control EcoRV's performance is surprisingly high. This 
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Figure 1: Sketch of an E.coli cell with native DNA (violet) 
concentrated in the centre. EcoRV restriction enzymes occur 
in two isomers: inactive, with closed cleft (red), and active with 
open cleft (yellow). Invading, foreign DNA (red double- helix) 
is attacked by EcoRV and cut. The void intracellular space 
shown here in reality is crowded by larger biopolymers. 



provides a concrete solution to the subdiffusion-efficiency 
paradox and supports current ideas that subdiffusion does 
not contradict efficient molecular reactions in cells. 

EcoRV enzymes and our model approach. 

The type II restriction endonuclease EcoRV binds non- 
specifically to double-stranded DNA. Once it locates 
its specific six-base sequence 5'-GAT|ATC-3' it cuts the 
double-strand and renders it inactive. This is an import- 
ant mechanism in the cellular defence against alien DNA 
stemming from, e.g., viruses attacking the cell [19]. The 
cell's native DNA is protected against EcoRV by mcthyla- 
tion of the DNA at cytosine or adenine . Interestingly, 
EcoRV is found in two configurations [2T][22]: as seen by 
X-ray crystallography, the unbound protein may switch 
between an inactive structure with a closed cleft and an- 
other, in which the cleft is more open. In the open, active 
state EcoRV binds both non-specifically to DNA and spe- 
cifically to its cognate binding site. 

In Fig. Q] we sketch an E.coli cell and its native DNA, 
EcoRV enzymes in active and inactive states being either 
attached non-specifically to the native DNA or freely dif- 
fusing in the cellular cytoplasm. The apparent void space 
in reality is a highly crowded ('supcrdcnse' [§]) complex li- 
quid, in which the enzymes subdiffuse. An invading DNA 
is being recognised an cleaved by active EcoRV enzymes. 

Remarkably the probability x act to find the enzyme in 
the open-cleft, active state at a given instant of time is 
as low as ~1% |23j . It is a priori puzzling why a vital 
defence mechanism should be equipped with such a low 
activity A physiologic rationale of the open/closed iso- 
merisation could be to reduce non-specific binding to the 
cell's native DNA. Alien DNA invading the cell would thus 
immediately be surrounded by a higher EcoRV conccntra- 




Figure 2: Diffusion trace (yellow) on a bond percolation cluster. 



tion that, after switching to the active state, could attack 
this DNA [23]. Our results show, however, under the as- 
sumption of normal diffusion in the cell the performance 
of EcoRV is only marginally better than that of a 100%- 
active mutant: normal diffusion on the length scales of an 
E.coli cell provides very efficient mixing, and the reduced 
activity of EcoRV would not constitute an advantage. As 
we will show this situation changes drastically under sub- 
diffusion, and low activity in fact becomes advantageous. 

Following previous studies [27H30] we model the subdif- 
fusion of EcoRV enzymes as random walks on a critical 
percolation cluster. The use of a fractal medium is in 
line with observations that the crowded cytoplasm may 
have a random fractal structure [25| . The lattice spacing 
is chosen as the size of an EcoRV. A bond between lat- 
tice points is created with occupation probability p. If 
p = 1 the entire lattice is accessible and the diffusion is 
normal. Reducing p to the critical value p c = 0.2488 or 
slightly above, the resulting cluster of permitted bonds is 
fractal with dimension df w 2.58 and the diffusion be- 
comes anomalous with a w 0.51. From extensive lattice 
simulations (see Appendix) we sample the times an en- 
zyme needs to locate its target, a specific sequence on an 
invading stretch of DNA randomly positioned in the cel- 
lular cytoplasm (the volume not occupied by the native 
DNA). The average target knockout time, equivalent to 
the mean first passage time (MFPT) to hit the target in 
an active state, is studied as function of the bond occu- 
pation probability p and the non-specific binding constant 
K° s of active EcoRV to DNA. 

Fig. [2] depicts part of the spanning percolation cluster 
and part of a random walk trace on this cluster. Due to the 
reduced connectivity in the fractal cluster, the trajectory 
reflects the existence of holes existing on all scales. The 
appearance of dead ends and bottlenecks in the scale-free 
environment effects subdiffusion [4]. 
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Figure 3: Typical time for the restriction enzyme to locate the 
target in an active state (MFPT) on a cubic lattice, as function 
of the bond occupation probability. Close to criticality (p c — 
0.2488 is marked by the vertical line), subdiffusion emerges 
with anomalous diffusion exponent a = 0.51. The non-specific 
binding constant is K% s = 10 7 [M -1 bp -1 ] [26]. Error bars are 
of the size of the symbols or less. 

Simulations results. — In Fig. [3] wc compare the 
MFPT for EcoRV (activity x act = 0.01) and mutant en- 
zyme (x act = 1) versus the bond occupation probability 
p, ranging from full occupation (p = 1, normal diffusion) 
down to the percolation threshold p — 0.2488 (subdiffu- 
sion with a = 0.51). For normal diffusion (p = 1) the 
MFPT is just a factor of two smaller for EcoRV, compared 
to the fully active mutant. Approaching the percolation 
threshold the native EcoRV increasingly outperforms the 
mutant, at criticality EcoRV's MFPT is two orders of mag- 
nitude shorter than that of the mutant. 

On average, the concentration of EcoRV is approxim- 
ately 1/iL'act = 100 times higher in the cytoplasm out- 
side the volume of the cell's native DNA than that of 
the mutant enzyme. At criticality, it is time-costly to 
cover distances, and thus EcoRV is 100 times more effi- 
cient than the fully active mutant enzymes. The latter 
become trapped around the native DNA, to which they 
bind non-specifically. In contrast, under normal diffusion 
conditions {p = 1) spatial separation is hardly significant, 
and the lower concentration is compensated by the higher 
activity of the mutant. 

In absolute numbers, even under severe anomalous dif- 
fusion with a w 0.51 EcoRV's MFPT is only a factor of 
ten higher than at normal diffusion. That means that the 
low-activity property of EcoRV renders their efficiency al- 
most independent of the diffusion conditions, compared 
to the huge difference observed for the mutant. The 
highly increased relative performance of the native EcoRV 
is the central result of this study. It demonstrates that 
subdiffusion is not prohibiting efficient molecular reac- 
tions. Moreover, our result provides a novel rationale for 




Figure 4: Sketch of the cross section of an E.coli model cell. 
Region 1 contains the cell's native DNA. In Region 2 ("cyto- 
plasm"), foreign target DNA are attacked by active restriction 
enzymes. The various symbols are explained in the text. 

EcoRV's elusive low-activity, that in this light appears as 
a designed property. We note that the MFPT shown here 
is the result for an individual EcoRV enzyme. Typically, a 
bacteria cell combines a fairly large number of restriction 
enzymes of various families. This significantly reduces the 
time scales indicated here, while preserving the character- 
istics of the EcoRV superiority. 

Discussion. — To obtain a better physical under- 
standing of the role played by events of EcoRV non-specific 
binding to the cell's native DNA, we study the depend- 
ence of the MFPT on the non-specific binding constant 
K® s . Experimentally, K® s can be varied by changing the 
salt concentration of the solution. For the cubic lattice 
(p = 1) we obtain an analytical expression for the MFPT 
for the geometry sketched in Fig.SJ We distinguish Region 
1 containing the native DNA, and Region 2 representing 
the cytoplasm, in which the foreign DNA enters and the 
EcoRV action occurs. 

Let us first address the non-specific binding of EcoRV 
enzymes to the native cellular DNA in Region 1, corres- 
ponding to a volume of length of L and radius R. As- 
suming rapid equilibrium with respect to enzyme binding 
and unbinding from the DNA, we observe the following 
relation between the volume concentrations of bound and 
unbound active (ready-to-bind) enzymes, 



In our notation c^ 1 ' is the overall volume concentration of 
enzymes in Region 1, while d and cj^ lk , respectively, 
measure the volume concentrations of enzymes bound to 
the native DNA and of unbound enzymes. The non- 
specific binding constant K® s to DNA refers to active 
(open-cleft) enzymes per DNA length, and is of dimension 
[K° s ] = M _1 bp _1 . Finally, Z^na is the length of DNA 
per volume in Region 1. Of the unbound enzymes, a frac- 
tion x act is in the active (open-cleft) state, ready to bind 
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to DNA. Thus, the concentration of active unbound en- 
zymes in Region 1 becomes c a ]) t = Xnctc^iki and one may 
introduce an overall binding constant K ns = x ac tK® s : 



c (1) 

L bulk 



1 



1 



^act-^Ss^DNA 



"ns^DNA 



(3) 



°bound 

As the total enzyme concentration in Region 1 is c- 1 ' = 
cl]L + cVL„a, and we have c ( ^ und = c^Zact-^L^NA we 
= cW/[l + - ^° 



Xacti^n^DNA]. In Region 1 



J bulk ' °bound' a 

c (1) 

the enzyme concentration in the continuum limit will be 
governed by a diffusion equation of the form 



dt 



(4) 



where D c ff = A3d/(1 + £act A^^na) i s an effective dif- 
fusion coefficient incorporating the assumption of rapid 
equilibrium with respect to binding to DNA and switch- 
ing between active and dormant states. ID diffusion along 
the DNA is assumed so slow that it can be ignored in con- 
nection with the overall diffusion of the enzyme. Indeed 
the ID diffusion constant for EcoRV have been measured 
to be orders of magnitude smaller than for 3D diffusion 

In Region 2 we assume that 3D diffusion is fast allowing 
us to write a conservation law for enzymes in the form of 
the difference between the flux across the boundary with 
Region 1, and the amount of enzymes reacting with the 
target per time, 



^ F (2) c (2) = _ A (l) Def[ 

at or 



k a c 



(2) 



(5) 



Here is the volume of Region 2, = 2ttRL is the 
surface area of Region 1, and k a is the rate constant for 
reaction with the target. The x act dependence of the rate 
k a is k a = ^act^a, where fc^ is the rate constant for the 
active state. Note that in our approach we assume that 
the switching between active and dormant state is fast 
in comparison with the diffusion across the regions, i.e., 
we may assume an equilibrium between these two states. 
Finally, we take Cbuik to be continuous across the boundary 
between the Regions 1 and 2, and that initially the system 
is at equilibrium with respect to the reaction-free situation 
with k a = 0. From the above system of equations the 
average search time yields in the form 



T = 1 



■iact-fins'rjNA 



a x act 



R 2 



1 



8Asd i + y 



(0) 



where y = V™ / [V^(l + x act K° s l^ 

The simulations were carried out on a 100 x 100 x 100 
cubic lattice with native DNA occupying a 100 x 50 x 50 
lattice in the middle. To compare the present calcula- 
tion with the simulations we assume a lattice spacing of 
a = 10 nm and set L = 1 fxm, = 1 fxm 3 and thus 
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Figure 5: MFPT on a normal lattice (p = 1), as function of 
the non-specific binding constant K® B . Simulations results are 
compared to the theoretical result ([6]). The vertical line marks 
the value K° B = 10 7 [M^bp -1 ] used in Fig. [3] 



= 0.25 ^m 3 . From this we obtain R = y/VW/(wL) « 
0.28 /im. Furthermore we choose the enzyme diffusiv- 
ity Z?3d = 3/im 2 /sec [value obtained for lac repressor 
in vivo at short times [35]] the DNA length per volume 



(i) 



= 1.5 x 10" 



l DNA 

that K„* = 



'm/yW, and K° s = lO^^bp -1 , such 
10 5 M _1 bp _1 when x act = 0.01 [33] (except 
for plots where this parameter is varied) with the base 
pair length bp = 0.35 nm. For the target association rate 
constant we take k a a = (Na) 3 /T lattico = a 3 (l - i? 3 d)/T s te P - 
Here T lattice = N 3 T stcp /(l - R 3d ) is the average search 
time for a random walker starting far from the target on a 
N x N x N cubic lattice and spending a time T s t ep per step 
to nearest neighbour sites. i?3d ~ 0.340537 is the walker's 
return probability to its origin [34]. Matching r s to P with 
the above diffusion constant through the mean squared 
displacement of the walker we obtain T st c P = a 2 / (GD^) ps 
5.6 fj,s. The assumption of a one lattice site target gives a 
target size of a = 10 nm. We have taken this target size 
to be lower than the in vitro effective sliding length [33] 
at optimal salt conditions, partly due to possible blocking 
on the DNA by other DNA-binding proteins. 

The above numbers yield J^4>NA ~ 3 x 10 5 , 
2/k ct =i « 10- 5 , j/k ct =o.oi « 10~ 3 , W/k° a » 2 sec, 
and i? 2 /(8£>3d) ~ 0.003 sec, and with these parameters 
we have to a good approximation T = K^J^^V^ /k°, 
regardless of whether x ac t = 1 or x ac t = 0.01. 

The simulations results for the case of normal diffusion 
are displayed in Fig. [5] At small A'° s the mutant (x ac t = 1) 
clearly outperforms EcoRV, the gap in the MFPT cor- 
responding to the reduced activity (x ac t = 0.01). At 
increasing K® s both EcoRV and mutant perform almost 
identically, with a small advantage to EcoRV. In this re- 
gime almost all active enzymes are bound to the cellu- 
lar DNA, such that EcoRV has approximately a factor of 
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l/^act = 100 higher bulk concentration. Concurrently, 
its association rate constant with the target DNA in the 
cytoplasm is reduced by the same factor. In this nor- 
mally diffusive regime dominated by non-specific binding, 
reduced activity of the restriction enzyme has no signific- 
ant advantage. The resulting MFPT behaviour according 

to Eq. © T w IOdna^^Ao m this regime depends 
linearly on the non-specific binding constant. Indeed, this 
behaviour is independent of x act . The agreement between 
the theoretical model and the simulations results is excel- 
lent over the entire range of K ° s (Fig. [5]). 

Fig. |6] shows the behaviour in the case of subdiffusion: 
almost over the entire K® s range, EcoRV significantly out- 
performs the mutant. At sufficiently large K® s values 
(above some 10 3 M _1 bp _1 ) the value of the MFPT is 
approximately two orders of magnitude smaller, i.e., the 
performance is improved by a factor close to the value 
1/a^act- This behaviour is thus dominated by the costly 
subdiffusion from the site of non-specific binding to the 
target. At low K® s values both curves converge. Now, 
the MFPT is fully dominated by anomalous diffusion to 
the target. Due to the compactness of the diffusion on 
the fractal cluster the difference between EcoRV and the 
mutant becomes marginal: upon an unsuccessful reaction 
attempt, EcoRV has a higher probability to hit the target 
repeatedly before full escape, improving the efficiency. In 
Fig- El the thick lines show the average of simulations over 
three different critical percolation clusters, while the thin 
black lines depict the result for each individual cluster. 
Apart from the low K® s limit, the results are very robust 
to the shape of the individual cluster. It would be interest- 
ing to derive analytically the MFPT dependence on K® s . 
While the MFPT problem on a fractal has been solved re- 
cently |35j , it is not clear how to apply this method in the 
present case, due to the division of the support into two 
subdomains. Similarly, for the related case of fractional 
Brownian motion |36j this remains an open question. 

Conclusions. — It is often argued that molecular pro- 
cesses in the cell could not be subdiffusive, as this would 
compromise the overall fitness of the cell due to the slow- 
ness of the response to external and internal perturbations. 
Here we demonstrated a solution to this subdiffusion- 
cfficicncy paradox: specific molecular design renders the 
efficiency of EcoRV enzymes almost independent on the 
exact diffusion conditions. Even though EcoRV are not 
always ready to bind, under subdiffusion conditions the 
low enzyme activity represents a superior strategy. 

Cellular subdiffusion may also be modelled by fractional 
Brownian motion (FBM) or continuous time random walks 
(CTRW) [37]. FBM shares many features with diffu- 
sion on fractal structures, e.g., the compactness and er- 
godicity. The essential observations for the MFPT found 
herein should therefore be similar for the case of FBM. 
In contrast, for CTRW subdiffusion the high probability 
of not moving in a given period of time will significantly 
enhance the advantage of EcoRV over the mutant: while 
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Figure 6: MFPT on a percolation cluster close to criticality 
(p c = 0.25) versus the binding constant K® a - The value K° s = 
10 7 M _1 bp _1 used in Fig. [3] is marked by the vertical line. 
Thick coloured lines: average over three different percolation 
clusters. Thin black lines: results for the individual clusters. 

being trapped next to the target EcoRV would have ample 
chance to convert to the active state and knock out the 
target. It will be interesting to compare our results to 
simulations on a dynamic percolation cluster. 

Subdiffusion-limitcd reactions generally increase the 
likelihood for biochemical reactions to occur when the re- 
actants are close- by [9l[T6]. Such a more local picture of 
cellular biomolecular reactions in fact ties in with the ob- 
served colocalisation of interacting genes [T7]. In higher 
cells the similar locality is effected by internal compart- 
mentalisation by membranes. It will be interesting to ob- 
tain more detailed information from single particle track- 
ing experiments in living cells, in order to develop an in- 
tegrated theory for cellular signalling and regulation under 
crowding conditions in living cells. 

* * * 
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Appendix. — Our simulations of the search process 
were carried out on a 100 x 100 x 100 cubic lattice. We 
considered bond percolation, i.e., for each pair of nearest 
neighbour sites a bond is constructed with a probability 
p. The searching random walker is only allowed to walk 
between connected sites (compare Fig. [2]). The largest 
cluster of connected sites was chosen and the remaining 
sites discarded. Of these remaining sites those within the 
central 100 x 50 x 50-size part constitute Region 1 with 
the native DNA. The number of such sites is denoted by 
and the volume of this region is thus VW = a 3 N^. 
The remaining sites outside this region constitute the cyto- 
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plasm. The number of these sites is N^ 2 \ and the corres- 
ponding volume becomes = a?N^ 2 \ A single target 
site is chosen randomly among the cytoplasm sites. 

The initial position of the searching random walker 
is chosen by first deciding whether it is bound or un- 
bound. With probability (V~W + V^)/[V^ + V^(l + 

^act^njfjNA)] ^ * s cnosei1 to be unbound. In this case its 
initial position is chosen randomly among all 
sites, and with probability x act it is chosen to be in the 
active state (otherwise it is in the inactive state). If the 
walker is initially chosen to be bound, then it is placed 
randomly among the native DNA sites, and its state is 
initially set as active. 

The initial time is set to t = 0, and the search is carried 
out according to the following algorithm: 

1. If the walker is in the inactive state or is situated in 
the cytoplasm, a time r step is added to the total time 
t. If the walker is active and situated in the region 
with the native DNA, then a random time is added 
to t. This latter random time is taken from an expo- 
nential distribution with average r stop (l + ^^4>na)- 

2. One of the 6 directions possible on a cubic lattice is 
chosen at random. If a bond exists to a neighbouring 
site in this direction, the walker moves to this site, 
otherwise it stays put. 

3. If the walker is in the active state, it switches with 
probability 1/5 to the inactive state. If instead the 
walker is in the inactive state, it switches to the active 
state with probability x a ct/5/(l — x ac t)- In the case 
when £ a ct = 1, this switching is turned off and the 
walker always stays in the active state. 

4. If the walker is on the target site and in the active 
state, the target is considered to be found, and the 
time t is recorded as the search time. Otherwise the 
iteration goes back to step 1. 

This procedure is repeated 5,000 times on the same per- 
colation cluster, but with a new target position and initial 
position of the searcher each time. The MFPT is calcu- 
lated as the average of the 5,000 recorded search times. 

References 

[1] See, for instance, P. H. von Hippel and O. G. Berg, J. 
Biol. Chem. 264, 675 (1989). For a recent example in- 
vestigating EcoRV, see Ref. [23] and M. A. Lomholt et 
al., Proc. Natl. Acad. Sci. USA 106, 8204 (2009). 

[2] A. P. Minton, J. Cell Science 199, 2863 (2006); H. Zhou, 
G. Rivas, and A. P. Minton, Ann. Rev. Biophys., 37, 375 
(2008); J. A. Dix and A.S. Verkman, Ann. Rev. Biophys., 
37, 247 (2008). 

[3] S. B. Zimmerman and S. O. Trach, J. Mol. Biol. 222, 599 
(1991). H. Dong, S. Qin, H.-X. Zhou, PLoS Comp. Biol. 
6, el000833 (2010). R. J. Ellis and A. P. Minton, Nature, 
425, 27 (2003). R. J. Ellis, Trends Biochem. Sci. 26, 597 



(2001). S. B. Zimmerman and A. P. Minton, Annu. Rev. 
Biophys. Biomol. Struct. 22, 27 (1993). 

[4] D. ben-Avraham and S. Havlin, Diffusion and Reactions 
in Fractals and Disordered Systems (Cambridge Univer- 
sity Press, Cambridge, UK, 2005). 

[5] C. D. Lorenz and R. M. Ziff, Phys. Rev. E 57, 230 (1998). 

[6] S. R. McGuffee and A. H. Elcock, PLoS Comput. Biol. 6, 
el000694 (2010). 

[7] G. Seisenberger et al., Science 294, 1929 (2001). 

[8] M. Weiss, M. Eisner, F. Kartberg, and T. Nilsson, Bio- 
phys. J. 87, 3518 (2004). 

[9] I. Golding and E. C. Cox, Phys. Rev. Lett. 96, 098102 

(2006) . 

[10] I. Bronstein et al., Phys. Rev. Lett. 103,018102 (2009). 
[11] S. C. Weber, A. J. Spakowitz, and J. A. Theriot, Phys. 

Rev. Lett. 104, 238102 (2010). 
[12] J.-H. Jeon et al., Phys. Rev. Lett. 106 048103 (2011). 
[13] D. Banks and C. Fradin, Biophys. J. 89, 2960 (2005). 
[14] W. Pan et al., Phys. Rev. Lett. 102, 058101 (2009). 
[15] J. Szymanski and M. Weiss, Phys. Rev. Lett. 103, 038102 

(2009). 

[16] G. Guigas and M. Weiss, Biophys. J. 94, 90 (2008). 

[17] G. Kolesov et al., Proc. Natl. Acad. Sci. USA 104, 13948 

(2007) . 

[18] M. A. Lomholt, I. M. Zaid, and R. Metzler, Phys. Rev. 

Lett. 98, 200603 (2007); I. M. Zaid, M. A. Lomholt, and 

R. Metzler, Biophys. J. 97, 710 (2009). 
[19] A. Jeltsch, C. Wenz, F. Stahl, and A. Pingoud, EMBO J. 

15, 5104 (1996). 
[20] A. Pingoud, M. Fuxreiter, V. Pingoud, and W. Wende, 

Cell. Mol. Life Sci, 62, 685 (2005). 
[21] S. G. Erskine, G. S. Baldwin, and S. E. Halford, Biochem. 

36, 7567 (1997). 
[22] F. K. Winkler et al., EMBO J. 12, 1781 (1993). 
[23] B. v. d. Broek et al., Proc. Natl. Acad. Sci. USA 105, 41, 

15738 (2008). 

[24] A. D'Arcy et al., J. Biol. Chem. 260, 1987 (1985). 

[25] R. Cuthbertson, W. M. L. Holcombe, and R. 

Paton, Computation in cellular and molecular biological 

systems (World Scientific, Singapore, 1995). 
[26] See Supplementary Material. 
[27] M. J. Saxton, J. Chem. Phys. 116, 203 (2002). 
[28] C. Loverdo, O. Benichou, and R. Voituriez, Phys. Rev. 

Lett. 102, 188101 (2009). 
[29] C. C. Fritsch and J. Langowski, J. Chem. Phys. 133, 

025101 (2010). 

[30] Similar insights into the effect of crowding environments 
were obtained from lattice random walks and off-lattice 
Brownian dynamics simulations in J. D. Schmit, E. Kam- 
ber, and J. Kondev, Phys. Rev. Lett. 102, 218302 (2009); 
N. Dorsaz et al., Phys. Rev. Lett. 105, 120601 (2010). 

[31] I. Bonnet et al., Nucleic Acids Res. 36, 4118 (2008). 

[32] J. Elf, G.-W. Li, X. S. Xie, Science 316, 1191 (2007). 

[33] A. Jeltsch and A. Pingoud, Biochemistry 37, 2160 (1998). 

[34] B. D. Hughes, Random Walks and Random Environ- 
ments, Vol. 1 (Oxford University Press, Oxford, UK, 
1995). 

[35] S. Condamin, O. Benichou, V. Tejedor, R. Voituriez, and 

J. Klafter, Nature 450, 77 (2007). 
[36] J.-H. Jeon, A. V. Chechkin, and R. Metzler, EPL 94, 

20008 (2011). 

[37] For a recent summary, see, e.g., S. Burov, J.-H. Jeon, R. 



p-6 



A solution to the subdiffusion-efBciency paradox 



Metzler, and E. Barkai, Phys. Chem. Chem. Phys. 13, 
1800 (2011). 



p-7 



